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TURBULENCE AND TEMPERATURE FLUCTUATIONS BEHIND A HEATED GRID 
By R. R. Mills, Jr., A. L. Kistler, V. O'Brien, and S. Corrsin 


SUMMARY 

In the approximately isotropic velocity and temperature fluctuation 
fields behind a hot grid, measurements were made of fluctuation levels 
and of various double and triple correlation functions. 

The double and triple correlation coefficient functions are of 
roughly the same spatial extent for the vector and scalar fields. As 
anticipated from theoretical considerations, the temperature fluctua- 
tions die out more slowly than does the turbulence. 


INTRODUCTION 

The simplest turbulent flow is that for which the statistical prop- 
erties of the field are invariant to rotation or reflection of the Car- 
tesian coordinate system. This notion, isotropic turbulence, was intro- 
duced by Taylor in 1935 (ref. 1) and has been fruitful in permitting 
relatively detailed analytical studies (ref. 2) . 

An approximation to isotropic turbulence turns up In the high Reyn- 
olds number flow far behind a plane "porous" obstacle spanning a uniform 
mean flow In a duct. The customary obstacle is a square mesh biplane 
grid of round rods. 

Although the turbulence found in natural and technological flows is 
ordinarily far from isotropic, many features of diverse turbulent flows, 
especially the spatially local features, seem to be moderately universal 
(ref. 3). 

A comparison between -the correlation equations for concomitant (in- 
compressible) isotropic velocity and temperature fields shows differences 
attributable to the fact that velocity Is a vector, while temperature is 
a scalar (refs. 4 to 6) . (As pointed out In refs. 4 and 5, the work 
applies equally well to isotropic turbulent mixing between two different 
constituents, provided the molecular mass transfer coefficient is 
nearly constant over the concentration range encountered.) The relative 
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decay rates were estimated under strong simplifying assumptions, and the 
relative "microscales" -were deduced in terms of the Prandtl number of 
the fluid. Since these theoretical predictions are essentially conjec- 
ture, it is necessary to determine some facts through measurement. 

A comparison of the corresponding statistical properties of these 
two fields has not only intrinsic interest; it may also contribute even- 
tually to an understanding of the difference between momentum and heat 
transport rates in turbulent shear flows with mean temperature gradients. 

A "hot grid" with thermal mesh equal to momentum mesh was selected 
with the expectation that approximate equality of integral scales would 
be obtained. The analysis of reference 4 suggests that the flow with 
equal integral scales for velocity and temperature fields may not be the 
simplest case analytically but is relatively simple for experimental 
realization. 

This work conducted at The Johns Hopkins University has been spon- 
sored and supported financially by the National Advisory Committee for 
Aeronautics. Acknowledgment is made to Dr. L. S. G. Kovasznay, Mr. L. ' 

T. Miller, and Mr. J. L. Luraley for their advice and Mr. S. Bhaduri for 
preparing the figures. 

SYMBOLS 

°p 

V*x> 

®u^*x) 
e,e^,eg 

f(r) 


g(r) 


i 


specific heat at constant pressure 

one-dimensional power spectrum of 0-fluctuations 

one -dimensional power spectrum of u- fluctuations 

output voltage of hot wire due to turbulent fluctuations; 
subscripts 1 and 2 distinguish between different wires 

longitudinal double vel ocity correlation coefficient, 

u(x,y,z)u(x-Kr,y,zJ _ u l u 2 
u’(x) u 1 (x+r) u^Ug 

l ateral double veloc ity correlation coefficient, 
ulx,y,z)u(x,y+r,zT 

u ,2 (x) 

wire current 

output of cubing circuit, microamps 
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K 

k(r) 


L f ,Lg 
M 


(r) 


m 


n(r) 

Pe 

Pr 

p(r) 

Re 

r 

S„ 


3 Ae 


S u 

S- 

u 

S Au( r ) 

S 6 


output of squaring circuit, microanqps 
thermal conductivity coefficient 
longitudinal triple velocity correlation. 


2, 
u ( 

[x,y,z> 

z 

lu(x+r,y,z) u l u 2 


[u- (xj] 

2 u’ (x+r) u’ 2 u^ 


integral scales of f and m, respectively 
grid mesh size for momentum and heat, M = 1 in. 


temperature correlation coefficient, -■ ^ ~ g^x+rf = ~§T~f 


Ci 0^ Q 

triple temperature correlation, — -= 

|_0'(x)_p 0' (x+r) 0^0^ 

Peclet number, u’Aq/V 
Prandtl number, v/r 

temperature-velocity triple correlation, 

u(x,y,z)0(x,y,z)0(x+r.y,zJ _ ^l 8 !^ 

( x ) 0'(x) 0 '(x+r) uj0£0£ 

Reynolds number , A u'/v 

s 

space interval 

skewness of voltage, e^/(e 2 )^ 2 


skewness of voltage derivative 


3/2 


skewness of 


’ fe) / p)_ 

voltage difference, (e-j_ - e 2^^{j- e 1 " e 2^ 


skewness of velocity, u^/(u 2 )^^ 


skewness of velocity derivative, (bu/bx)^ j J^(3u/Sx)^j 
skewness of velocity difference, (uj_ - Ug ) 3 /Jj u i " u 2 ) 2 ] 3/2 


skewness of temperature fluctuation, ft 3 / ( 6 2 ) 


3/2 
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n 

T mean fluid temperature _ 

t time 

U mean flow velocity — 

u,v turbulent velocity fluctuations in x- and y-directions , 

respectively 

u-l u(x) 

u 2 u(x+r) 

V 1 = v ( y )> v 2 = v ( y+r ) 

x,y,z Cartesian space coordinates , x alined with mean flow 
X thermal diffusivity coefficient, K/pCp 

T|,T|g length parameters in universal equilibrium theory for velocity 

and temperature, respectively (t) is the "Kolmogoroff micro- 
scale” ) 

© temperature fluctuation _ 

©^,©2 = ©(x) , ©g = ©(x+r) 

x c ,x 0 wave numbers at lower and upper bounds of inertial subrange, 
respectively, corresponding to velocity spectrum except 
when subscript 0 is used 

x x wave number in x-direction _ 

7\g,'A^.,Ag dissipation scales and lateral and longitudinal velocity and 
* ^ temperature fluctuations, respectively 

ji viscosity coefficient 

v kinematic viscosity coefficient, p/p 

£ space interval in x-direction, = Ax 

p density 

Subscripts : 

max maximum 


V 


wave number corresponding to temperature spectrum 
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Superscript : 

indicates root -mean- square values , e. g. , u' = 

EXPERIMENTAL EQUIPMENT 

Approximately isotropic turbulence was produced at seme distance 
behind a biplane square mesh grid of l/ 4-inch-round metal rods spaced 1 
inch at the centers. The grid was positioned in the wind tunnel as 
shown in figure 1. The signals generated by the hot-wires in the wind 
tunnel were fed into the circuit outlined in figure 2. In order to 
measure the various skewness factors and triple correlation functions 
presented later in this report, it was necessary to incorporate the 
power and operational amplifier circuit (fig. 3) to boost the signal 
levels to the proper input level of the squaring and cubing circuits 
shown in figures 4 and 5. (The auxiliary squaring and cubing circuit 
was developed by Mr. L. T. Miller of the Aeronautics Department, The 
Johns Hopkins University. ) 

A typical set of calibration curves for this auxiliary squaring and 
cubing circuit is shown in figure 6. The response of the auxiliary am- 
plifiers was flat to approximately 10 kilocycles when used as ordinary 
power boosters j and when used as differentiating operational amplifiers, 
the response was linear with frequency to approximately 6 kilocycles 
(fig. 7). 

The wind tunnel walls were adjusted to provide uniform mean flow 
over approximately 90 percent of the cross- section. This configuration 
produces a turbulence which is not quite isotropic: v’ 0.9 u' at 

x/M = 100. (This is consistent with all previous measurements in which 
the authors have participated (see, e.g., ref. 7). 

The grid was heated with 220-volt, 3-phase alternating current in 
order to generate approximately isotropic temperature fluctuations in 
this flow field. No 120-cycle-per-second periodic temperature component 
was detectable with this arrangement. To produce a uniform mean temper- 
ature field downstream of the grid, it was necessary to have local con- 
trol over the heating current supplied to grid rods. Adjustable resis- 
tors were added to the grid heating current supply lines to achieve the 
desired mean temperature distribution (fig. 8) . 

The present arrangements generated satisfactorily homogeneous tur- 
bulence and temperature fields (fig. 9). 

Turbulence and temperature root-mean- square fluctuation measurements 
were made with the hot-wire anemometry equipment described in reference 
8. The experimental procedures for s imultaneous velocity and temperature 
fluctuations are given in reference 9. 
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All the measurements were made at a mean speed of 14 feet per sec- 
ond and a mean temperature rise (across the grid) of about 5° C. The 
resulting temperature fluctuations were inconveniently small, but this 
ensured negligible influence of density variations upon the fully devel- 
oped turbulent velocity field, as attested by both turbulence level and 
velocity correlation measurements. 

Because of the extremely small temperature fluctuations, it was nec- 
essary to use unusually high resistance hot-wires, 0.00005-inch plati- 
num of about 75- to 150-ohm resistance and about 0.1 inch long. No 
length corrections were applied to the data. The high resistance neces- 
sitated correction for heating current fluctuations when the system was 
operated as an anemometer. The uncompensated time constant of these 
wires was about 0.1 millisecond. In this particular flow field, the 
wire current giving equal root-mean- square signal contributions from ve- 
locity and from the temperature fluctuations was 5.5 milliamperes. The 
sensitivity to temperature was approximately linear with current, and 
the sensitivity to velocity was proportional to the current cubed for_ 
currents less than roughly 0.7 milliampere. 


EXPERIMENTAL PROCEDURE 

The decays of the temperature and velocity fluctuations were meas- 
ured with a single hot-wire used at various operating conditions so that 
the separate decay functions could be calculated (ref. 9). The correla- 
tion was zero with the accuracy of the data. 

Two-point double correlations were measured with two wires mounted 
on a traversing device that located the wires with respect to each other 
within 0.005 inch. One wire was held stationary and the other was moved 
•downstream, and the line joining the wires made an angle of 5° with the 
mean flow velocity. This was done in order to avoid interference effects. 
When the wires were at their closest position, they were laterally sepa- 
rated rather than one behind the other. Appropriate corrections have 
been applied to the correlation data by assuming the isotropic relation 
for the correlation tensor. 

Since preliminary measurements showed that the velocity double cor- 
relations were unaffected by the presence of tenperature fluctuations, 
the velocity measurements reported here were done with the heat off for 
convenience . 

Temperature double correlations were measured using the hot-wires 
as resistance thermometers (i^ 1 ma). The correlations were corrected 
"for noise by assuming the noise uncorrelated with the signal. 
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For convenience, the u’ spectrum, was measured without heating 
after it had been determined that the u' levels and correlations were 
not appreciably changed by the heating. 

The temperature spectrum was measured with "mixed" sensitivity, 
that is, wire temperature set for the same order of response to velocity 
and temperature. This was done because, under pure-resistance-thermometer 
operation, the signal was inconveniently low at the higher frequencies. 

Of course, increasing the wire temperature does not increase the rate of 
gathering of information on the larger voltage signal. 

With mixed sensitivity the hot-wire responds to both velocity and 
temperature spectra, which must be separated. Simple superposition of 
the two energy spectra follows from the assumption that all harmonics 
of the two spectra are uncorrelated. Since 0u 0 in this field, such 
an assumption seems reasonable. Then the temperature spectrum is ob- 
tained by subtracti n g the (unheated grid) velocity spectrum from the 
mixed spectrum. 

The skewness of the velocity spatial derivative was determined by 
differentiating the velocity signal with respect to time and measuring 
the mean cube and mean square of the resulting signal. The proportion- 
ality of instantaneous space and time derivatives is close for this flow 
(refs. 10 and 11). Since the noise was uncorrelated with the derivative 
signal, the noise power could be subtracted directly from the mean-square 
signal. The noise skewness was measured and found equal to zero so that 
no correction was made to the cubed signal. 

The triple velocity correlation measured in this experiment was 
k(r) since k(r) requires hot-wires sensitive only to u velocity. The 
wire locations are, of course, identical to those required to measure 
P(r). 

To measure u^Ug(i) the signals from two hot-wires, separated a 

distance | = r, were subtracted and cubed. The resulting signal, when 
assuming identical hot-wires, is proportional to 

3 3 , _ 2 , 2 

u-l - u 2 + ou^Ug - ou-^Ug 

If the sum of the two signals is also cubed, the result Is 

3 3 2 2 

Uf + Ug + 3u-^Ug + 2 > u 1 u 2 

For isotropic turbulence, u^Ug ( § ) is an odd function of | so that 


u f u 2 = " 
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Also 


u i = 4 = 0 

Therefore , if isotropy is assumed, the difference signal is just 6u^u 2 , 

the desired quantity, and the sum is identically zero. When the sum was 
measured for the wires in their closest position, however, it was found 
that the mean cube of the signal was eight times the mean cube of a 
single-wire voltage. This anomalous behavior was traced to the nonlinear 
distortion of the signal by the hot-wire itself. The analysis of this 
situation is given in the appendix. 


When the difference of the two voltages is used, the nonlinear hot- 
wire effect on the mean cube is considerably smaller, a phenomenon sim- 
ilar to the canceling of second-harmonic distortion by a push-pull am- 
plifier. Therefore, all measurements of the correlation u^Ug reported 


here are measured by using the difference signal only. 


The function s AuW * defined as (u-j_ - U 2 ) 


| = r, was obtained by direct measurement of the 
both numerator and denominator go to zero but in 
ratio is equal to the skewness of the derivative 




>) 


3/2 


for 


s a1°> * 

signal. 


u 2 , 

two factors. At r = 0, 
such a way that their 
(ref. 2). Therefore, 


was obtained hy direct differentiation of the velocity 


The cube of the difference signal from two wires adjusted to be 
equally sensitive to velocity and temperature fluctuations is given by 


(e^ - eg)^ = u^ - ug + Sugu^ - 3u^ug + 0^ - ©g + 30§0q ~ 3®l% + 

^ U 1®1 + 3^1 “ ® U 1®1®2 “ 3ug0^ - 3ug0g + 6ug0j_0g + 

p O - p p ■ ■ ■ ■ — 

39^u£ + 30-^Ug - 60-^u-^Ug - 30gU^ - 3©2 U 2 + (l) 


With isotropy, all the terms containing the variables at only one point 
and the velocity to an odd power are zero. There are also the following 
isotropic relations : 
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U 2 U 1 = " u l u 2 


*l e 2 = U 2 0 l 


^Vl = "iVb 

Hie correlation u^Sg or Ugfi^ is identically zero for an iso- 
tropic, incompressible field just as the pressure velocity correlation 
is zero (ref. 10 ). The correlation u-j_u 2 ©]_ is equal to u^u 2 0 2 since 

this is one component of a second-order isotropic tensor, solenoidal in 
one index, just as f and g are even in the pure velocity case. 
Therefore, continuity gives a relation between the two components, and 
this relation preserves the symmetry property. The components analogous 
to g (e.g., S-^v-^Vg) are obviously symmetric, and, therefore, Q-^u^Ug is 
also symmetric in £. 

With the assumption of isotropy in the field, therefore, the aver- 
age of the difference voltage cubed is finally reduced to 


u i e i S 2 ^ 2^2 


e, e r 


e l e 2 - Vi 


fl 3 - A 3 
e l - e 2 


(e-j_ - eg)^ = 6u^u^ + 12u^0-^0 2 (2) 

To determine the triple temperature correlation 0£0 2 ( f=), the two 

wires must be sensitive only to 0 temperature fluctuations. In this 
measurement, both the sum and difference signals of the two wires must 
be cubed. By assuming identical hot-wires, the cube of the sum minus 
the cube of the difference signal is proportional to 


0^0 2 (ep + e 2 ) 5 - (e 2 _ - e 2 ) 3 (3) 

which is the required correlation. 

The preceding analysis assumes identical hot-wires. In practice, 
this is an extremely difficult condition to meet. It was found much 
easier to determine the sensitivities of the two wires by actual cali- 
bration and to modify the above analysis to account for this diff erence. 

2 

The modification is straightforward in the determination of u-^ug and 
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0^© 2 5 however, in the measurement of the calculations require » 

that each wire he equally sensitive to velocity and temperature fluctua- 
tions. This was impossible to accomplish in general, and the measure- 
ments of this function contain contributions from correlations neglected 
by symmetry conditions in the preceding analysis. However, rough esti- 
mates of the order of magnitudes of these terms were made, and appro- 
priate corrections have been made. 


EXEERIMEnmL RESULTS 

The velocity and temperature fluctuation decay curves are plotted 
in figures 10 and 11, respectively. A typical pair of runs is given for 
each to indicate the repeatability of the data as well as the equality 
of hot- and cold-grid flows. The inverse squares are plotted in figures 
12 and 13. 

The time spectra of 9 and u at x/M =17.0 were determined as 
described in the previous section and are shown in figure 14. These may 
be considered fair approximations to the one-dimensional longitudinal 
space spectra, except for the lowest wave number range (ref. 11 ). With- 
in the experimental precision there is no significant difference. The * 

circuit noise has been subtracted out for both sets of data. 

The skewness of the velocity derivative, defined by 



is shown in figure 15 as a function of x/M. For isotropic turbulence, 

S. = k'" (0)7v^. The curve indicates that it is a slowly increasing func- 

tion of x with values in the vicinity of 0.4. The nonconstancy may 
indicate that the Reynolds number of the turbulence was not sufficiently 
high for to be determined solely by eddies in the Kolmogoroff in- 

ertial subrange (ref. 3) ; there were contributions from the larger eddies 
and from those influenced by viscosity. Other comparable data, including 
some measurements by Batchelor and Townsend at the same Reynolds number 
behind a grid of lower solidity (ref. 12) and one measurement by Stewart 
(ref. 13) , are included for comparison. 

The temperature field was produced by the turbulence-producing grid 
itself. It was hoped that this method of introducing the temperature 
fluctuations would not cause a ’’ spotty" condition to be generated. 
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However, figure 16 shows S Q = ET / (d ) as a function of x/m. Al- 
though S g is small (max. value, ~ 0.08), this evidence of an initial 
"pulse" character is quite measurable. 

The longitudinal velocity correlations are shown for various values 
of x/ M in figure 17. The definition of this correlation is 
uiu 2 

f(r) = — ; — - where the subscript 1 denotes spatial location x, y, z and 
^1^ 

subscript 2 denotes x+r,y,z. Primes on u-j_ and u 2 mean root-mean- 
square value. It is assumed that f(r) is essentially equal to the iso- 
tropic correlations of von Karman and Howarth (ref. 10). If the turbu- 
lence were similar in structure at all decay times, a plot of this f(r) 
against a nondimens ional scale proportional to a local characteristic 
length would collapse all curves into one. It is well established that 
such is not the case in actual turbulence. For this plot the length 
selected was the "dissipation scale" or "Taylor microscale," defined by 



The A’s used here were obtained by fitting parabolas to the vertices. 
Independent dete rm ination was made via the decay equation valid for iso- 
tropic turbulence. This equation is also due to Taylor: 


.si = -I0v 1 

at >1 

and with the use of Taylor's hypothesis for the interchangeability of 
space and time on the average: 


du 2 v u 2 v u 2 

= -10 \ = -20 % 

^ u h 2 u 

g f 


( 5 ) 


Therefore, A can be deduced from the decay curve for u 2 as a function 
s 

of x. The comparison is given in figure 18, along with one value 

obtained from the second moment of the measured power spectrum (ref. 2) : 




( 6 ) 
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Within the experimental scatter, the use of r/Aj* collapses the 

correlations in the region near the origin hut not at the tails, consis- 
tent with earlier experimental results. The systematic decrease of the 
correlation with decay time at a fixed r/A f shows the well known fact 

that the turbulence is not accurately similar; still, it is ro ugh ly so. 
The temperature double correlations ;~ 

8-t So 

' = 9 . 

were plotted (fig. 19) in the same way as the velocity correlations but 
against r/7 Vg, where this thermal "microscale" is given by (ref. 4) 



That the temperature fluctuation field was reasonably isotropic is 
evidenced by the approximate equality of the longitudinal and lateral 
temperature double correlations (fig. 19(b)). 

The isotropic temperature decay equation gives an independent esti- 
mate of Ag (ref. 4) : 

= -is x. £, 
dx U a| 

The Ag scales are compared in figure 18 along with those (ref. 5) com- 
puted from the temperature spectrum by an equation like (6) . 

Figure 20 shows a typical contrast between f (r) and m(r) at the 
same station. Although Ag ~ A^, which corresponds to rough equality 

of vertex curvatures, m(r) runs higher immediately thereafter, which in- 
dicates that 

m iv( 0 ) > fiv(o) (9) 

The experimental curves for the triple velocity correlations are 
shown in figure 21. The values of k(r) were always negative for the” 
x/M and r values covered. According to isotropic turbulence theory, 
the triple correlation is zero at the origin, rises with an initial 
cubic* variation, and is zero at r = °°, The signal fluctuation made it 
impractical to measure at large r where k » 0. Whether k(r) changes 


( 8 ) 
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* 


* 


sign or goes to zero with a negative power greater than -4 (ref. 2) can- 
not he determined from these data. 

Figure 22 shows mixed temperature-velocity correlation coefficient 
at two of the x/M values where the triple velocity correlations were 
measured. As measured. 


p(r) 


u l®l e 2 



( 10 ) 


These data are somewhat more scattered than those of k(r) . When a cor- 
relation function was measured several times, there sometimes occurred 
systematic differences attributed partly to the inevitable difference 
between any two hot-wires. Furthermore, individual wire properties 
change as the wire is used, mainly because of the accumulation of dirt 
from the airstream. 


As shown previously, under the assumption of isotropy it is possible 
to get p(r) without knowing anything except k(r) if the wires are 
identical and have their velocity and temperature sensitivities equal. 

In practice, however, if the inevitably nonidentical wires are operated 
so that they have the same sensitivity to velocity, their temperature 
sensitivities will be different. This condition necessitates measurement 

o 

of other correlations (such as or< 3- er 1° obtain p(r). Esti- 

mates were made of these corrections to p(r) an d wer e negligible for 

the three stations presented. (The correlation 8^9 Z was available from 
a separate measurement. ) 

The temperature-velocity correlation coefficients have a shape sim- 
ilar to that of the velocity correlations and are of the same order of 
magnitude. They were negative for all x/M and r studied. 

Figures 21 and 22 show the triple correlation behavior as a func- 
tion of decay time. These curves are plotted against an abscissa scale 
selected to collapse the curves in the neighborhood of the origin. Anal- 
ysis, if the fields are considered as isotropic, shows that the velocity 
correlation starts in the well known way: 


k(r) 


u x_ 
3 3.' 
g 


(H) 
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and that 


p(r) 



( 12 ) 


Therefore., since is roughly constant, k(r) plotted against r/A-p 

should he universal near the origin. Whether p(r) should be a universal 

function of r(7\ f 7\|) for r -*■ 0 is determined by whether 


■ 8 f is roughly a constant as is the skewness factor of 

u' 0 2 X 2^ “ 2 

Bu/ax. The plots indicate that f \"^) 


increased more rapidly 


u' 0“ ^ ^ . 2 

with decay time than S^. The quantity l"5cj vas no ^ m sasured di- 
rectly because of the loss of the signal in the tube noise. 


It appears from figure 23 that p__„ increased with increasing 

x/M for the Reynolds and Peclet numbers of this experiment. The peak 
values of k(r) also increase in absolute value for increasing x/M in 
agreement with the result of Stewart (ref. 13) . Near the origin they 
fall into one curve, a necessary consequence of the relative constancy 
of S^. This is evidenced in figure 21 (a) where equation (ll) has been 

computed and plotted. 

Plotted in this fashion the tails of neither p nor k are uni- 
versal functions, which again shows that the fields do not remain similar 
at different decay times. 

Figure 24 shows the function 

( U 1 ~ 

at several x/M. This function attains its maximum at y = 0 where its 
value is S^, the skewness factor of the velocity derivative (ref. 2) . 

It falls off rapidly to about 0.2 and. then gradually decreases, presum- 
ably to zero. Its behavior, at all x/M covered, is the same. The 
fact that this function is not a constant iB another indication that the 
Reynolds number of this turbulence was not sufficiently high for a dom- 
inant inertial subrange. These measurements agree very well with those 
of Stewart (ref. 13). 
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In figure 25 n(r) is shown as a function of r for the three sta- 
tions examined. Although the scatter In the data is greater in these 
results than in the other triple correlations presented, the general 
trend of this correlation is indicated clearly. The required limit, 
n(0) = Sg, is substantiated within the accuracy of these experiments. 

Figure 26 shows the measured values of g(r). Included for compar- 
ison is the function g(r) computed from the well known isotropic rela- 
tion (ref. 15) 

gCr) = h 'k [} 2f ^ 

The disagreement shown is probably partly due to the inaccuracy asso- 
ciated with any graphical differentiation process, partly because the 
scale characterizing the turbulence changes for r > 0 in the f (r) 
curve, and, fl.uaJ.ly, because the turbulence is not precisely isotropic. 


ANALYSIS 


The theories of isotropic turbulence and of isotropic scalar fluc- 
tuations in isotropic turbulence are still unsolved; in fact, they have 
not even been formulated in a determinate way without recourse to ad hoc 
postulates. Nevertheless, a number of consequences of isotropy have 
been inferred, including some particular (and relatively simple) forms 
of the averaged differential equations. Several of these consequences 
have been invoked in the preceding section. 


Two-point correlation equations are of special interest here. For 
isotropic fields decaying in time, there are for the vector and scalar 
fields, respectively (refs. 4 and 14): 


and 


since 






+ 2vu 2 


(£t 

[hr 2 


+ 


4 Bf\ 
r "Sr j 



^cu 2 ) 1/2 





k(r) = -2h(r) 
h(r) = - k(r) 


(14) 


( 15 ) 
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Sufficient data were taken to permit computation of every term in 
these equations. Since the triple correlations were the least certain 
experimental results , these measured functions were compared with those 
obtained from the double correlations by equations (14) and (15) . 

Figure 27(b) shows a typical comparison for the velocity field. 

For small r/?Vj> the agreement is somewhat less than that of Stewart's 
data (ref. 13) , which are plotted in figure 27(a), but not significantly 
so. Figure 22(b) gives the comparison for the mixed triple correlation 
at station x/M = 32.0. 

The small disagreement is presumably attributable to either lack 
of isotropy or inaccuracy in the measurement of triple correlations. 
Since the small structure should be reasonably isotropic even if the 
large structure is not, the latter reason seems more likely. 


If the correlations appearing in the above equations are expanded 
in a power series around r = 0 and the coefficients of the tP terms 
are kept, the so-called "decay equations" result: 



d6 2 

dt 


-I2r -2= 


(16) 


(17) 


The actual experiments are done in steady state behind a grid with 
the various averaged quantities depending upon x instead of t. Taylor 
has pointed out, in effect, that, if ' these "vary sufficiently slowly with 
x, _equation (16) can be applied to the grid flow with the substitution 
x/U = t. The same holds for equation (17), so that 



-20 -zz~ — p 

U Af 


ae 2 

dx 


-12 £_ if 

s 


(18) 


(19) 


These two equations permit calculations of A^(x) and Ag(x) from 

the experimental data on u' (x) and 8' (x) , as mentioned in the previous 
section. The faired results are included in figure 18. 
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Also plotted in figure 18 are the values of A obtained graphi- 
cally from the double correlation curves according to 


_1_ m” (o) 
7,2 2 
A e 

i_ _ t" Co) 

" 2 


( 20 ) 

( 21 ) 


and those calculated frcan the power spectra. 

Reference 1 shows that for two special "hinds" of isotropic fields 
the ratio ~Xq/7\£ has the same simple value: 



( 22 ) 


where Pr is the Prandtl number v/r- The two kinds of isotropic fields 
are: 


(1) Reynolds and Peclet numbers are so smell that the convective 
effects are negligible for both heat and momentum. Therefore, this is 
not really turbulence in the ordinary sense. 

(2) Reynolds and Peclet numbers are both very large, and each field 
is assumed to have completely "self -preserving" correlation function dur- 
ing decay. In this model proposed by von Karman (ref. 15), the dissipa- 
tion and conduction terms of the two correlation equations are neglected, 
but decay is included indirectly by substitution of the decay equations 
into the terms of the correlation equation. 


Equation (22) can also be deduced by an approach somewhat different 
from the previous: 


( 3 ) The Reynolds and Peclet numbers are assumed to be so high that 
the inertial subranges (with ® u ( x x ) ancl Eg(x x ) both are suf- 
ficiently extensive to permit the approximations 




« X 


• 5/3 


X 


x 0 <x x < 


= 0 


elsewhere 


( 23 ) 
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E e (x x ) 


** x 


= 0 


-5/3 


_ c Xy < x . 

0,© * c,0 

elsewhere 


(24) 


The cutoff wave number x c is just the inverse of Kolmogorof f 1 s 
microscale: 


2 * 

x c = T 


and x q is the analogous thermal quantity (ref. 5) : 


2k 


In reference 5 it is shown that 

x 


c ’ 6 V e 


= (Pr) 3 / 4 


(25) 


(26) 


(27) 


The microscales can be obtained from the second moment of the power 
spectrum as 




x J E u( x x) dx x 


x JC E e (x x )dx x 


(6) 


(28) 


Substituting equations (23) to (27) 
Xq^ 3 relative to x^/ 3 and Xq/q 


into (6) 
relative 


and (28) 

x c{§ 


and neglecting 
result in 



(29) 


If consideration is now restricted to fields In which the thermal large 
eddy structure is roughly equal to the velocity large eddy structure 

(as in these experiments), x n *» x n a) which again results in equation 

( 22 ). 0 0,6 
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In fact, none of the three sets of assumptions leading to equation 
(22) is directly applicable to these experiments. The Reynolds and 
Peclet numbers are not low enough for case (l) and not high enough for 
the other two cases. Furthermore, it was shown that the correlation 
functions of velocity do not remain completely similar during decay. 

The Important conjecture to he made is that a result encountered under 
such divergent sets of assumptions may have approximate validity in a 
wide range of situations, perhaps including the actual one. 

By taking Pr = 0. 72 for air, equation (22) gives 

^ « 1. 2 (30) 

A f 

which is in reasonable agreement with the downstream values of Ag/A^ 

computed from the decay curves but is higher than the values computed 
from the correlations. The data indicate a trend toward satisfaction 
of this relation for large decay times. The only other measurement of 
this ratio was on the axis of a heated round turbulent jet (ref. 8), 
where the value was 0.95. 


The solid Ag/A^ curve in figure 18 is just a particular way of 

displaying the decay data. Two other ways are also instructive. Figure 
28 shows the variation in Reynolds and Peclet numbers : 


Re s 
Pe = 



(31) 


with A taken from the decay data. The apparent constancy of Pe is 
evident. 


Figure 29 shows the relative decay rates of temperature and veloc- 
ity fluctuations. For small x/M, these fluctuations die out at about 
the same rate, but, for larger x/M, 9' dies out at a decidedly lower 
rate than does u 1 . For fully isotropic -fields the ratio of equation 
(19) to equation (18) gives 



du 1 ” 10 (Pr) \A,J 
u' 


( 32 ) 
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and, If equation (22) holds, 


ML _h1 _ 3 

6 ' du' ~ 5 


(33) 


which is very close to the measured results at large x/M. Notice that 
equation (33) is independent of Pr. 


No attempt has been made to compare the individual decay rates of 
9 ' and u-' with the theoretical predictions that can be obtained by a 
variety of postulates. 


The integral scales are defined by 


and 



f (r)drl 



m(r)dr 




(34) 


These are shown in figure 30. Reference 4. shows that, with self- 
preserving correlations, if the scales are equal at any instant and also 
at the initial time, they will remain equal. The small discrepancy here 
is doubtless a result of the fact that there is certainly no complete 
similarity in either field. Another consequence of a supposed similar- 
ity in both fields is that 


(35) 


Figure 31 shows good agreement with both relations for large decay times. 

Figure 23 shows the value of k^^ as a function of decay time, 
contrasted with the results of Stewart (ref. 12) . It is interesting to 
speculate on the history of for even greater decay times. A mono- 

tonic increase in k (toward unity) would not be in contradiction to 

max 

the usual approximations for very large decay time , that is , the so- 
called "final period" results, because the coefficient of the k term of 
the correlation equation is Hie factor which certainly becomes negligible 
as Re -*■ 0. It seems plausible, though, that in the limiting case of 


and that 


7 — Pe = const. 
h 9 


'A 


> 


~ Re = const 
L f 
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random Stokes' flow, analogous to heat conduction in a solid medium., the 
velocity probability density is symmetrical, in which case k(r) would 
go to zero at t -» “ . 


CONCLUDING REMARKS 

In roughly isotropic concomitant velocity and temperature fluctua- 
tion fields it has been found that: 

1. The temperature fluctuations die out more slowly than the turbu- 
lence (fig. 29) in approximate agreement with the predictions of refer- 
ence 4. 


2. The temperature double correlation coefficient function is rough- 
ly equal to the "longitudinal" velocity correlation. 

3. The "mixed" triple correlation coefficient function entering the 
double temperature correlation equation is of the same order as the cor- 
responding velocity triple correlation. 

More explicit statement of detailed quantitative results covered 
herein will not be given here because these results may not have suffi- 
cient universality. They may be significantly characteristic of the de- 
viations from isotropy in these velocity and temperature fields. 


The Johns Hopkins University, 

Baltimore, Md. , June 12, 1957. 



22 


NACA TN 4288 


APPENDIX - EFFECT OF HOT-WIRE NONLINEARITY 

ON TRIPLE CORRELATION MEASUREMENT 

At "small" fluctuation levels, perhaps u'/U < 10 percent, the 
Dryden-Kuethe linearized equation for hot-wire response (ref. 16) is 
satisfactory for measurements of a number of turbulence properties. For 
example, it is well known that turbulence level, double correlation, and 
power spectrum measurements are not seriously distorted in the (<3- 
percent level) isotropic turbulence rather far behind a grid. However, 
previous research (.ref. 13) makes no mention of this possible source of 
error in triple correlation measurements. It is evident that this re- 
placement of the response curve by its local tangent will lead to poorer 
approximations for functions which depend more heavily on the tails of 
the probability density of u(t) . Furthermore, away from response curve 
inflection points it seems evident that odd moments will tend to fare 
worse under linear approximation than will even ones. 

The simplest case of concern here is the skewness of u, 
s „ = » 5 /c u^ )3/2. A first-order est im ate of the error in this measure- 
ment is obtained by keeping the parabolic terms in the perturbation form 
of King's equation for hot-wire sensitivity. A suitable approximate 
form of King’s equation is 

i i 2 ^ p — 

= A - i 2 + bJu (Al) 

8 * 

where 

A,B constants 

i wire current, taken constant here 

R wire resistance 

R„ wire resistance at ambient fluid temperature, taken constant here 
U instantaneous velocity normal to wire 

Writing U = U + u and R = R + r where u = r = 0 and expanding the 

2 

square root out to u give 
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Averaging equation (A2) and subtracting the average from (A2) 
result in 



from which could be computed the measured skewness. 



as a function of the true skewness S u and of other flow and wire prop- 
erties. For simplicity, the analysis is restricted to isotropic turbu- 
lence, in which case S u = 0. The relation u 4 /( u 2 ) 2 = -3 is taken as 
the Gaussian value, and the final estimate is 


S 


e 




A - 



i 2 + B 



_u_ 

U 


(A4) 


For a typical 0.1-mil plat inum wire at an overheat temperature of 50° C 
in 2-percent turbulence level flow, this estimate gives 


** 0.05 
e 

even though S u = 0. 

The seriousness of this effect can be gauged by observing that 0.05 
is the order of the maximum value of the triple correlation coefficient 
k(r) . 


Before proceeding to an estimate of the effect on triple correla- 
tion measurement, the determination of skewness of velocity derivative 
is considered: 
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Following a proce dure similar to that for S e , an estimate is obtained 



For the typical hot-wire operation cited previously, equation (A5) gives 

S. « -S. (l - 0.005) 
e u v . 

which Indicates a negligible effect. The quadruple correlation is esti- 
mated from the data of reference 17. 

In a completely formal measurement of triple correlation coefficient 
(using mean cubes of sum and difference voltage) an additive error can be 
expected of the order of that in S e , which is appreciable. However, for 

isotropic turbulence the triple correlation is just proportional to the 
skewness of the velocity difference (see section on experimental proce- 
dure) , and it is not difficult to show that in velocity difference skew- 
ness measurement the error due to nonlinearity is small . In a sense, the 
effects of the two wires tend to cancel. Hence, it is convenient for 
reasons of both speed and reduction of error to follow Townsend in assum- 
ing isotropy and thereby using the ideal relation 

W “ u 2^ 3 = “ 6 u 1 u 2 = ~ 6 (u 2 ) ^ k (r) 

Writing the two wire responses in series form gives 

CD 00 

e i = y] a i u i e 2 = X) (as) 

i=0 j=0 

where and a^ are proportionality constants. Terminated with the 

u 2 terms, equations (A6) become equation (A3), from which an estimate of 
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the first-order error committed by assuming that any dimensionless sta- 
tistical function of e-j_,eg is equal to the corresponding function of 


u^Ug. In particular, the difference skewness factors 

( e l - eg^/jje-L - eg) ^ and (u x - Ug) 3 /^^ - u 2 ) 2 J ^ 
to be related by 


S A = S A <1 - 3 
As Au 


*(•♦*) = 


o n 4 ~2"1> 

1\ ^ “ u l u 2 


U 2 (u-j_ - u 2 )‘ 


|(o 2 + o + |)-^- 


7t2 /■ \ 2 

U K - u 2 } 


turn out 


higher order terms 


(A7) 


For the example used earlier equation (A7) gives an error of less 

than 0.2 percent for r > ^ A and rises to perhaps 2 percent at r -*■ 0. 

Hence, this effect is negligible within the accuracy of these 
measurements . 


In the cubed- difference approximation to triple correlation, the 
same approach leads to 



after the assumption of U'/Cu 2 ) 2 = 3. This corresponds to percent errors 
about 10 times those in S Ae' but these errors are appreciable only as 


r 0 where (e-^ - eg 
no difficulty arises 


l — Mia ii 

\3/ „2 . 2 

' / 1^/ 2 


is also going rapidly to zero; thus. 
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Figure 2, - Block diagram of electronic components. 
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Figure 4. - Squaring circuit. (Ail tubes 12AT7.) 
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Figure 6. - Concluded. Typical calibrations of squaring and 
cubing circuits. 
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(a) Auxiliary and hot-wire amplifiers. 

Figure 7. - Frequency response of measuring equipment. 
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Figure 7. - Concluded. Frequency response of measuring equipment. 
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Figure 9. - Homogeneous turbulence and temperature fields. 
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Time spectra of u 2 and 0 2 . x/jM = 17.0. 


Figure 14. - Concluded. 
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(b) Actual data measured at x/M « 17.0. 

Figure 19. - Concluded. Temperature double correlations for various x/M, 
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(b) x/n * 32.0. 

Figure 21. - Continued. Triple velocity correl&tiona for various x/k. 
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Figure 25. - Triple temperature correlation n(r) as function of r. ' ’ 
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Figure 26. - Measured 
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(a) Data of reference 12. x/M = 60.0. 
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(b) Present data. x/M. = 32.0. 

Figure 27. - Comparison of measured values of triple correlation with 
values predicted from equation (14). 
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